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[ Abstract. We introduce hermitian and unitary quantum phase operators of two- 

component fermionic fields. The phase operators of bosonic fields are the same as that 
^ , of quantised radiation fields introduced long back. Using these phase operators we 

study quantum phase and number fluctuations of a pair of interacting bosons and those 



of a pair of interacting two-component fermions in a symmetric double-well potential. 
i <~| The fluctuations in inter-well number- difference operator can be quantified in terms 

of two-mode squeezing parameter which also serves as an entanglement parameter to 
quantify inter-well entanglement in number variables. We investigate the effects of on- 
site interaction on the fluctuations of the two canonically conjugate inter-well number- 

■ and phase-difference variables which are closely related to the entanglement. We also 
, study intra-well number fluctuations. Our analytical and numerical results reveal that 

the fluctuation properties of a pair of bosons are quite distinct from that of a pair of 
two-component fermions. These findings may have implications in possible quantum 
information processing with ultracold atoms in confined geometries such as double- well 

■ optical lattices. 
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1. Introduction 

Quantum phase fluctuation of interacting systems plays an important role in describing a 
variety of physical phenomena such as phase transitions, superfluid tunneling, Josephson 
effects etc. With advancement in precision interferometry with ultracold atoms [TJ |2j [3] 
in confined geometries such as traps and atom chips, developing a proper understanding 
of the quantum phase properties of interacting many-particle systems is of prime 
interest. As far as quantum phase measurement and its theoretical interpretations 
are concerned, there are unresolved issues which need to be addressed, in particular 
in the context of emerging field of atom optics. For instance, a proper definition 
of 'quantum phase' of electromagnetic fields had remained a hotly debated topic in 
theoretical quantum physics for a long time |31 El El IH] ■ Accurate determination of 
phase difference between two optical fields in the quantum domain remains an elusive 
task due to lack of theoretical understanding of quantum phases. It is also known that 
such a difficulty exists even in the case of semi-classical radiation theory if the field is 
weak; and amplitude and phase fluctuations are correlated. About two decades ago, 
Mandel's group [91 [TUl [HI [12] experimentally examined two closely related but distinct 
measurement schemes for determining phase difference between two optical fields in 
both semi-classical and quantum cases. They made use of sine and cosine of phase 
difference operators as defined by Carruthers and Nieto [7] . Furthermore, to resolve the 
problems associated with quantum phase measurement, Noh et al. [91 [11] introduced 
an operational definition of quantum phase that requires different phase operators for 
different measurement schemes. 

Quantum phase problems for electromagnetic fields had been extensively studied 
during 90's [13]. However, in the context of matter waves or atom optics, quantum 
phase has not been addressed adequately. Physical situations in case of matter waves of 
ultracold atoms are different from that in electromagnetic radiation. Trapped ultracold 
atoms can be considered as an isolated interacting many-particle quantum system in 
the absence of any appreciable trap-loss. It is then necessary to formulate quantum 
phase of interacting matter waves with fixed total number of particles. This makes a 
major difference when one compares quantum phase of light with that of matter waves 
since the photons of light are noninteracting and their total number is not fixed. To 
understand the underling features of quantum phase of an interacting massive system, we 
resort to a simple microscopic model of two particles interacting via a contact potential 
in the ground state of a symmetric 1-D double-well. There are several advantages 
of this model. First, this enables us to work with an analytical solution. Second, 
considering the fact that only phase difference between two modes of an electromagnetic 
field is measurable, analogously the ground states of two wells of a double-well can be 
regarded as two modes of matter waves. Third, this model allows us to study the effects 
of exchange symmetry of the two particles on quantum phase properties in a simple 
analytically tractable manner. This model can be readily generalised for optical lattice 
where one can study inter-site quantum phase fluctuations and their effects on Bose- 
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Hubbard physics [21 [02 [TBI HZ] and superfluid-Mott quantum phase transition [16J. 
Furthermore, since symmetric double-well potentials are important for revealing the 
importance of symmetry in a wide range of quantum dynamical systems starting from 
rotations of molecules to spontaneous symmetry breaking in quantum field theory, this 
model of quantum phase fluctuations may be connected to many areas of research. Single 
particle dynamics using ultracold atoms in a double-well has been studied in different 
contexts such as quantum tunneling [HI [19], Josephson oscillations [18] , nonlinear self 
trapping [201 EI] of Bose-Einstein condensates, coherent control of tunneling [22] etc. 

In recent times, quantum dynamics of an interacting atom-pair in a double-well 
has become important. Because, several experiments with ultracold atoms in a double- 
well potential have revealed significance of interaction in correlated atom-pair tunneling 
[21 [23], number squeezing [3] in a Bose Einstein condensate (BEC) and exploration of 
quantum magnetism [21] with cold atoms. Theoretical studies with interacting atoms 
in a double-well have demonstrated entanglement in atomic hyperfine spin and phase 
variables [2H], an interplay between interaction and disorder in a BEC [2E], operation 
of a quantum gate [2Z] and so on. With increasing use of double-wells in cold atom 
research, double-well optical lattice [211 I2H1 |3Ql [31] is emerging as an important tool 
for studying correlation effects in cold atoms [32J. With this perspective in interacting 
atoms in double- wells, quantum phase fluctuations are important for gaining new insight 
into correlation dynamics of cold atoms in such a confined geometry. 

In this work, we address quantum phase problem for interacting bosons and a pair 
of interacting two-component fermions. In addressing this problem, we realise that, 
although the definition of a hermitian and unitary quantum phase of massive bosons 
will remain the same as that of photons, a proper quantum phase operator for fermions 
seems to be not available in the literature. Since in reality, only phase differences 
are measurable and not the absolute phase, we first introduce hermitian and unitary 
quantum-mechanical phase- difference operators of two-component fermions. In our 
double- well model, total number of particles is a constant of motion. In order to elucidate 
canonically conjugate nature of number-difference and phase-difference operators, we 
apply our formalism to study quantum fluctuations of both number- and phase- difference 
operators. One can introduce two non-commuting operators corresponding to the cosine 
and sine of the phase- difference operators. Both of them are canonically conjugate to 
the number- difference operators. These two phase operators plus the number-difference 
operator form a closed algebra. A unitary and hermitian phase-difference operator 
and corresponding phase-difference state can be constructed using the cosine and sine 
phase-difference operators. We find that the fluctuation in cosine of phase-difference 
of a pair of two-component fermions is greater than that of two bosons for otherwise 
same parameters for all times and for any arbitrary on-site interaction strength U. 
For both the bosonic and fermionic cases, the fluctuations in cosine of phase-difference 
tend to maximize with diminishing number-difference fluctuation in the limit U — > oo 
for times t « h/U where h stands for Planck's constant. We also study on-site 
fluctuations to understand the bunching or anti-bunching properties of two particles 
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in a double-well and their possible connection with quantum phase fluctuations. On- 
site number fluctuation of two non-interacting bosons shows periodic anti-bunching 
behavior. In contrast, two non-interacting two-component fermions do not exhibit any 
number fluctuation. As interaction is switched on, in the low interaction regime, on- 
site fermionic number fluctuations show intriguing oscillatory bunching property with 
modulations, collapses and revivals. Collapses occur at some specific times depending 
on the ratio of U to the tunneling strength J. Counter-intuitively, on-site number 
fluctuation of an interacting pair of two-component fermions never exhibits any anti- 
bunching. In contrast, on-site number fluctuation of two interacting bosons shows 
strong anti-bunching behavior in the regime of low U. Oscillatory bunching behavior 
of two-component fermion-pair in the low or intermediate interaction regime implies 
interaction-induced correlated pair tunneling and consistent with similar tunneling effect 
for two-component bosons as experimentally observed by Bloch's group [21 [23]. The 
dynamical equations for a pair of two-component bosons and those of a pair of two- 
component fermions are the same. The reduced phase fluctuation in two-component 
fermions or two-component bosons as compared to that in single-component or spin- 
polarized bosons might be related to the oscillatory bunching or correlated tunneling in 
two- component case. 

This paper is organised in the following way. In section 2, we present the formalism 
of two-mode quantum phase- difference operators. The model and its solution are 
described in section 3. The results are discussed in section 4. We conclude this paper 
in section 5. 

2. Quantum Phase Operators 

Before we discuss our model problems, we present here operator formalism that is 
particularly suitable for describing quantum phase fluctuations of bosons and fermions in 
a double-well. For bosons, a proper quantum mechanical phase operator can be defined 
following that for quantised radiation fields. For fermions, there exists no standard 
definition of a proper quantum mechanical phase operator. We here introduce a phase 
operator for fermions - this operator is hermitian and unitary. Before we discuss our 
new formalism, let us digress little bit and have a revisit into the history of quantum 
phase problem. 

Dirac |4J first postulated the existence of a hermitian phase operator in his 
description of quantised electromagnetic fields. However, Dirac's phase operator is not 
unitary. Susskind-Glogower [6] introduced a new hermitian phase operator which was 
again non-unitary. Louisell [5j first introduced the periodic operator function in defining 
a phase variable conjugate to the angular momentum. Carruthers and Nieto [7] showed 
that one can define two hermitian phase operators C and S corresponding to cosine and 
sine of the classical phase, respectively. Using these operators, they introduced two- 
mode phase difference operators of a two-mode radiation field. Explicitly, the two-mode 
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Figure 1. Behavior of the fluctuation ACf 2 of the cosine of phase-difference (left 
panel) and the fluctuation AVF B of number-difference (right panel) as a function of 
dimensionless time t — Jt for different values on-site dimensionless interaction strength 
U = U I J. From top to bottom, the 4 subplots in both the panels correspond to U = 
(a,e), U = 0.5 (b,f), U = 20 (c,g) and U = 500 (d,h), respectively. In the limit U -> oo, 
for t « 1/U we have ACf 2 ~ ^ an d AW B = meaning that number fluctuation 
is minimised while phase fluctuation is maximised in the large interaction regime in 
conformity with number-phase uncertainty. 



phase- difference operators are defined as 

C" CN /nf I c C 

12 — < -^l<-'2 T 0\02 

S™ = SxC 2 - S 2 C X (1) 

where 

Ci = i((JYi + iyhi + aj(iv 4 + (2) 

Si = ^(iV, + 1)-^ - aliN, + (3) 

are the phase operators corresponding to sine and cosine, respectively, of the i-th mode. 
In terms of creation (annihilation) operator a\(a,i) of the corresponding modes (z = 1, 2), 
phase difference operator can be written as 

C?2 N = + l)-hxal(N 2 + + a\(Nx + l)^(N 2 + l)~h 2 ] 

S?2 N = ^[{Ni + l)-h ia l(N 2 + - + + 1)-U 2 ] (4) 

The above cosine and sine phase-difference operators are non-hermitian and non- 
commutating. Pegg and Burnett [8] first introduced a hermitian and unitary phase 
operator. By synthesizing the methods of Pegg-Burnett [8] and Carruthers-Nieto, Deb 
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et al. [33] introduced hermitian and unitary phase-difference operators of a two-mode 
field with fixed number of total photons. What they did was to couple the vacuum state 
of one mode with the highest Fock state of the other to introduce the phase- difference 
operators of the form 



where 
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describe the contributions from the vacuum states of the two modes. \N\,N — N\) 
represents a two-mode Fock state with N\ and N — Ni being the photon numbers in 
mode 1 and 2, respectively. In case of quantised electromagnetic fields, the assumption of 
a fixed number of photons is made to circumvent the problem of non-unit arity. However, 
after all calculations are done one has to take the limit that the number of photons goes 
to infinity. 



2.1. Quantum phase operators of bosons in a double-well 

Here we discuss quantum phase operators for massive particles in a double-well potential. 
Since it is possible to keep the total number of particles in a double-well fixed in the 
absence of any loss, the assumption of a fixed total number of quanta (in this case, 
the total number of particles) is a good one unlike that in electromagnetic fields. We 
consider that the two modes 1 and 2 correspond to the second quantised matter wave 
fields to the left and right well, respectively, of a double-well potential, and the operators 
ctj(aj) represents annihilation (creation) operator of a particle in left (i = 1 and right well 
{i = 2). Then the operators defined in equations (5) and (6) suffice to be the two- mode 
phase-difference operators of bosons in a double-well potential corresponding to the 
cosine and sine, respectively, of the two-mode phase-difference. Thus we have Cf 2 = C12 
and Sf 2 = 5i2 where the superscript B stands for boson. The difference of the number 
operators or the population imbalance between the two wells is W B = a\a\ — a^a^. The 
three operators Cf 2 , Sf 2 and W B operators obey closed cyclic commutation algebra as 
follows 



[C? 2 ,S*),W l 



+ 



[Sf 2 ,W B ],Cf 2 + [W B ,Cf 2 ],S 



'12 







(9) 



The commutation algebra of the given operators are following 



Cf 2 ,W l 

sf 2 ,w B 



2i[S B -(N + l)S 
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Figure 2. Same as in figure 1 but for a pair of two-component fermions. From top 
to bottom, the subplots correspond to U = (a,e), U — 0.5 (b,f), U = 20 (c,g) and 
U = 500 (d,h), respectively. In the limit U -> oo, ACf 2 ~ ^ and AW F = 0. 



One can define a unitary phase-difference operator 

p 12 = Cf 2 + iS? 2 . (11) 

The eigenstate of this operator, that is, phase- difference state can be constructed as 
the product of single-mode phase states of Pegg-Burnett [8] subject to the condition 
that total number of quanta in the two modes is a constant of motion. The procedure 
for deriving phase- difference state is described in references [33j [M]. The important 
point to be noted here is that for low number of total bosons N, the effect of vacuum 
states such as | 0, N) and | N, 0) on quantum phase-difference would be significant. 
Therefore, in case of low N, one has to use unitary quantum phase-difference operators 
as defined in equation ffTT]) for accurate measurement of quantum phases. The phase- 
difference operators of Carruthers and Nieto will approach unitarity in the limit iV — > oo. 
In a recent theoretical paper by Sarma and Zhou [26], an operator similar to that of 
Carrruthers and Nieto has been implicitly used for studying phase dynamics of a Bose- 
Einstein condensate (BEC) in a double-well. It is worthwhile to mention that while in 
case of BEC, probably a non-unitary phase operator such as used in [26] can suffice for 
phase measurement for all practical purpose, a few bosons in a double well necessarily 
require unitary phase difference operators for high precision phase measurement, in 
particular for the purpose quantum information processing with a few bosons in a 
double-well. 
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2.2. Quantum phase operators of a pair of two-component fermions in a double-well 



Here we define quantum phase operator for fermions. The operators corresponding to 
the cosine and sin of two-mode phase-difference of a system of two-component fermions 
with spin up and down components can be defined as 

1 
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Here a^a^) is the fermion annihilation (creation) operator of spin component a in the 
z-th mode. Similarly, the other two operators Cf, 2 -t an d S'ii2t can ^ e defined. 
The fermion population imbalance operator is given by 

w F =J2 (fti* - ^) = E - a L« 2 .) (15) 

o-=t,4 cr=t,4 



The operators C*f 2 , S F 2 and H /F obey the closed commutation algebra 



+ 



[s! 2 ,w F ],c! 2 + [w F ,c! 2 ],s! 2 







(16) 



A unitary phase-difference operator for a pair of two-component fermions in the 
ground state of symmetric double-well potential can be written as 



12 



;i7) 



and corresponding eigen phase-difference state can be obtained in a similar manner as 
in bosonic case with N = 2. 



3. Model and its solution 



As a model system, we consider a pair of interacting particles in the ground state 
of a symmetric double-well potential with a finite barrier but infinite side walls. The 
interaction is assumed to be of zero-ranged contact type. We describe quantum dynamics 
of two cases: (1) A pair of indistinguishable bosons, and (2) A pair of two-component 
fermions. The bosonic model we use here is similar to the one studied by Longhi [SB] 
who has shown that the average dynamical behavior of a pair of hard-core bosons in a 
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Figure 3. The left panel shows the variation of fluctuation of the sine of phase- 
difference ASf 2 a s a function r for the bosonic case while the right panel exhibits the 
fermionic counterpart of the same. From top to bottom, subplots correspond to U = 
(a,e), U = 0.5 (b,f), U = 20 (c,g) and U = 500 (d,h), respectively. 



symmetric double-well has a classical counterpart in the transport of electromagnetic 
waves through wave-guide arrays. However, such an analogy can not be drawn for 
fermions in general. 



3.1. Two Bosons 

The Hamiltonian of a system of N interacting bosons occupying two weakly coupled 
lowest states of a symmetric double well potential is given by 



H = —KJ{a\a r + a\ai) + ^-(dj 2 d 2 + a} 2 a 2 r ) 



(18) 



where ai^{ai, r \) are the bosonic particle annihilation(creation) operators for the two sites 
/ (left) and r (right) of the double well, J > accounts for the hopping or tunneling 
between the two sites and U is the on-site interaction. The wave function \ip{t)) in the 
basis of Fock states with fixed total particle number iV can be written as 



N 



(19) 



^2 c ni {t)\n h N - m) 

ni=0 

Where c ni (t) is the probability amplitude to find n/ particles in the left well and (N — ni) 
particles in the right well and \ni, N — ni) denote Fock state with n\ particle at the left 
site and N — ni particle at the right site. From schrodinger equation , we obtain 

.dci(t) 



dt 



-(nici+i + + Vici 



(20) 
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where Ki = J\J (I + 1)(N — I), Vi = ^[l 2 + (N — I) 2 — N] and the normalization condition 

N 

is£|q(t)| 2 = l. 
1=0 

In the special case of N = 2, the solution is simple and analytically tractable. For 
instance, let us consider the initial condition c (0) = 1 and Ci(0) = c 2 (0) = i.e. both 
particles are initially in the right well. One then obtains 

Ci(t) = q(0) + y/2 

Cl (t) =Cexp [in B )r] +Dexp [-zOjr] (22) 

where subscript i stands for either or 2, U = U/ J, t = Jt, f2 B = \/4 + (U /2) 2 , 
Sl± = n B ± U/2, C = ^[fi B d(0) + v^MO) + Co (0)}], and D = ^[ft+ Cl (0) - 
V / 2{c2(0) + co(0)}]. Assuming initial condition co(0) = 1 and ci(0) = C2(0) = 0, i.e. 
both bosons are initially in the right well, we get C = ^=^- and L> = — ^-^ . Note that 
these solutions are the same as in [55] . 



C 



exp[if2 B r] — 1 ^exp[— iQ B r\ — 1 



Q1 



(21) 



3.1?. Two Fermions 

The dynamical behavior of two fermions in a symmetric double-well potential is 
described here. Since fermions obey Pauli exclusion principle, more than one particle 
can not be accommodated in a single quantum state. Here we deal only two-component 
fermions, one with spin up and another with spin down. The Hamiltonian of the system 
is given by 

H — — (a^a^a^a^ + a^ar^Q-^a^) (23) 

Now, we consider our trial wave function as a linear superposition of the Fock states as 
follows, 

W)) = ^(| t4) + I l,t» + c 2 (t)\ U,0) + c 3 (t)\0,U) (24) 

where the states | j, l) , | ^, f) define one fermion in the left well and another in right well 

and | f>lo 0) , |0,tl) define both fermions are in the left well and right well respectively. 

Here c\(t), c 2 (£) and c$(t) are the probability amplitudes of finding one fermion in one 

well, both fermions in the left and in the right well respectively. Now, putting equation 

(25) and (26) into Schrodinger equation, we get 

dci r- . U 

l ~dt = ~ ^ 2J(c 2 + c 3 ) + — ci 

i^=~ ^ Jc i + ^ ( 25 ) 

.dc 3 r- U 

%—r- = — v2Jci H c 3 

alt 2 3 

By solving equation (27), we get 

d(t) = Aexp [ifijr] + £exp [-i(fip r ] ( 26 ) 
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Cj (t) = Cj (0) + V2 



exp[z^ F r] - 1 _ ^expf-zfi^r] 



(27) 



where subscript j stands for either 2 or 3, Qp = y4 + (U /4) 2 , where Qp = Qp ± {7/4, 
A = ^[^(0) + v/2{c 2 (0) + c 3 (0)}] and B = ^[O+c^O) - v/2{c 2 (0) + c 3 (0)}]. 
Assuming initial condition 03(0) = 1 and Ci(0) = 02(0) = 0, i.e. both fermions are 
initially in the right well, we get A = ^ Q and B = — ^ n ■ 

4. Results and discussions 

Corresponding to the cosine and sine of the two-mode phase-difference operators for 
two bosons in the double-well, equation (10) governs the number-phase uncertainties 
between the number-difference and either of the two phase operators related to cosine 
and sine of the two-well phase- difference. Number-difference fluctuation AW B 2 is closely 
related to the the two-mode squeezing parameter which in turn describes entanglement 
between two bosonic modes in terms of number variables. Therefore, quantum phase 
fluctuation is also related to the two-mode entanglement. The entanglement between 
the two wells in number variables can be quantified as the two-mode squeezing [35] or 
entanglement [36], |37] parameter. 

£ n = (A(n, - n r )) 2 /((ni) + (n r )) = V N ) (28) 

The two wells (two modes) become entangled when £ n is less than unity, or 
equivalently, when becomes less than y/~N. For two bosons in the double- well, 

this parameter becomes 

2 . „,„ . 1 
I 

We display the temporal evolution of ACf 2 for different U values in the left panel of 
figure 1. We assume that both bosons are initially in the right well. The right panel 
of figure 1 shows variation of AI^ 6 as a function of r . The two wells are entangled 
if AW B < y/2. From fi gure 1 we notice that when bosons are noninteracting ({7 = 0), 
the number-difference fluctuation periodically oscillates between zero and v2 meaning 
that the two wells remain entangled for all the times. When Ajy B = 0, two wells are 
perfectly entangled. In noninteracting case, perfect entanglement varies periodically 
with a period of tt/2. With the increase of interaction strength U, the time periodicity 
of perfect entanglement increases. In the large U limit, we have AW B ~ for the times 
t « 1/U meaning that the two bosons become localised for long times. Now, in which 
well the two bosons will be localised depends on the initial condition. For instance, if 
the two bosons are initially in the right well, then the bosons will remain localised in 
the right well in the large U regime. As U increases ACf 2 start rapidly oscillating. The 
mean value of ACf 2 increases with the increase of U. In the limit U — > 00, ACf 2 settles 
down at a value y/3/2 for times r « 1/U. 



rli 
So 



1 ~ TF sin2 ^ BT ) ~ 4fi^ B cos ( fi B r ) - n B cos(fi B r)) 2 



(29) 
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Let us now discuss number and phase fluctuations in fermionic case. The 
entanglement parameter in one spin (a) component of fermions can be expressed as 

= 1 " 47? ^ C ° S( ^ T) " Q * COS ^)] 2 (30) 

The fermion entanglement parameter ^ in terms of total fermion number fluctuation 
is given by 



e F = 2 



2 1 

sin 2 (fi F r) - ^2"[^f cos(fijtr) - fi+ cos(fi F r) 

F F 



(31) 



The fluctuations of the phase- difference and population imbalance are given by 
(ACf 2 ) 2 =\ + \{c\ + c 2 ) + ^(c* 2 c 3 + c*c 2 ) - i( C ;c 3 + c*c 2 ) 2 

(AS F 12) 2 = \ + i(c 2 + c 2 ) - \(c*c 3 + c*c 2 ) + l -{clc, - c*c 2 ) 2 

(AW7 F ) 2 = 4 [(c 2 + c 2 ) - (c 2 - c 2 ) 2 ] (32) 

Figure 2 is a fermionic counter-part of figure 1 showing the phase-difference and number- 
difference fluctuations for the same parameters as in figure 1. While number-difference 
fluctuation shows almost similar temporal behavior as in bosonic case, phase- 

difference fluctuation ACf 2 is strikingly different from that of bosons. ACf 2 is always 
smaller than that of bosons for a particular value of U. However, in the limit U — > oo, 
AC*^ settles down at y/3/2 as in bosonic case. In the low interaction regimes, temporal 
evolution of number- difference fluctuations in both bosonic and fermionic cases shows 
intriguing interference pattern with modulations due to two times scales which are J 
and U. 

We plot the fluctuation of the sine of phase-difference AS*^ for both bosonic (q = 
B) and fermionic (q = F) cases in figure 3. We have verified that AC^ 2 and AS^ 2 satisfy 
the uncertainty conditions for both bosonic and fermionic cases. 

Next, in order to understand how on-site number fluctuations behave in bosonic 
and fermionic cases, we quantify the on-site number fluctuations in terms of Mandel Q 
parameter which is well known in quantum optics and serves as an important figure of 
merit for characterizing quantum statistical properties of light. For bosons in a double- 
well we define Q parameter as = ((An s ) 2 ) — (n s ), where s — l,r stands for either left 
(/) or right (r) well. Q characterises deviation of particle distribution from Poissonian 
statistics for which Q — 0. Sub-Poissonian (super-Poissonian) statistics corresponding 
to Q < (Q > 0) implies anti-bunching (bunching) of the particles. Let us Assume that 
both the bosons are initially in the right well. Then, in a very low interaction regime 
(U 0), we have Qf(r) ~ -| - 1 cos(4r) + cos(2r) and Qf(r) ~ -|sin 2 (2r). This 
implies that in the low interaction regimes, on-site Mandel Q parameter is always less 
than or equal to zero. Q < implies non-classical behavior. For 2r = (2n + l)n or 
equivalently r = (n + l/2)7r we have Qf ~ —2 and ~ 0. This means that at these 
specific times, number fluctuation in one well exhibits maximum non-classical behavior 
while that in the other well number fluctuation lies at coherent level (Q — 0). In a very 
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Figure 4. Mandel Q parameter Qi for on-site number fluctuation in the left well 
as function of r for bosons (left panel) and fermions (right panel) for U = (a,e), 
U = 0.05 (b,f), U = 0.1 (c,g) and U = 50 (d,h). 



high interaction regime (U — > oo), we find Qf ~ 1 — 16 ^ {/2 - > 0, — — 1 — 16 ^ [/ 2 - >■ — 2 
This means that the two bosons become localised in the right well for ever. In case of 
double-well lattice model, this will correspond to a Mott insulator phase with 2 bosons 
at each right lattice site. 

Number fluctuation of fermions in the s-th (s = l,r) well may be characterized 
by Qf( T ) = (tM'M) ~ (^stX^W)- m ver y ^ ow ~~ ^ 0) interaction regime, the 
fluctuation quantities behave as Qf (r) — >■ 0, (r) — >■ In very high ([/ — >• oo) 
interaction regime, all fluctuation quantities also go to zero. These analytical results 
are further corroborated by numerical results shown in figure 4 where we plot the 
variation Mandel Q parameter as a function of time for both fermionic and bosonic cases. 
Figure 4 shows that two-component fermions never shows any on-site anti-bunching. 
With increase of U, a pair of two-component fermions exhibit on-site bunching. In 
contrast noninteracting bosons shows periodic anti-bunching and coherent behavior, as 
U increases bosons exhibit both anti-bunching and bunching behavior at different times. 
Figure 4 indicate that while phase coherence (low U regime) in case of fermions may be 
related to coherent behavior, that in case of bosons appears to be related to non-classical 
non-Poissonian (anti-bunching) behavior. 

5. Conclusions 

In conclusion we have introduced quantum phase operators for two-component fermion 
fields in ground state of a symmetric double-well. Using the quantum phase operators 
we have studied in detail the effects of on-site interaction on quantum phase and number 
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fluctuation properties of two interacting bosons and a pair two-component fermions in a 
double-well potential. Our results reveal the effects of on-site interaction on dynamical 
evolution of squeezing in phase or number variables of a pair of interacting particles 
in a symmetric double-well system. This may be of interest to characterize quantum 
entanglement in continuous variables in confined geometries. A pair of two-component 
fermions interacting via a contact-potential in a symmetric double-well has average 
dynamical character that is qualitatively similar to that of a pair of bosons. Our study 
on two-mode quantum phase and number fluctuations reveal that these fluctuation 
parameters can clearly distinguish between two bosons from a pair of two-component 
fermions in otherwise similar physical conditions. In terms of number variables, both 
bosonic and fermionic systems exhibit interesting inter-well entanglement properties 
which may be a potential resource for future quantum information processing with 
neutral atoms in double-well optical lattices. Our results indicate that the correlated 
tunneling of a pair of atoms with two spin components as experimentally observed 
by Chen et al. [23] are related to bunching behavior with reduced inter-well quantum 
phase fluctuations. With the first demonstration of homo-dyne detection of a fluctuating 
continuous variable of matter waves by Gross et al. [39] last year, it might be possible 
in near future to perform experiments on the measurement of quantum phases of matter 
waves in a similar manner as in Mandel's experiments that require homo-dyne or hetero- 
dyne detection of weak signals. 
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